Frequency generation by a magnetic vortex-antivortex dipole in spin-polarized current 
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A vortex-antivortex (VA) dipole may be generated due to a spin-polarized current flowing through 
a nano-aperture in a magnetic element. We study the vortex dipole dynamics using the Landau- 
Lifshitz equation in the presence of an in-plane applied magnetic field and a Slonczewski spin-torque 
term with in-plane polarization. We establish that the vortex dipole is set in steady state rotational 
motion. The frequency of rotation is due to two independent forces: the interaction between the two 
vortices and the external magnetic field. The nonzero skyrmion number of the dipole is responsible 
for both forces giving rise to rotational dynamics. The spin-torque acts to stabilize the vortex dipole 
motion at a definite vortex-antivortex separation distance. We give analytical and numerical results 
for the angular frequency of rotation and VA dipole features as functions of the parameters. 



I. INTRODUCTION 

The dynamics of the micro-magnetic structure in fer- 
romagnetic elements presents unusual features (TJ [5]. 
Switching between different magnetic states which is a 
fundamental issue for magnetic recording, and magnetic 
oscillations which are exploited for the design of fre- 
quency generators [3] are examples for areas of appli- 
cations where magnetization dynamics is of primary im- 
portance. The most promising method to probe the dy- 
namics of magnetic structures is the injection of spin- 
polarized current in magnetic nanoelements. 

Of particular interest, in this respect, are magnetic vor- 
tices which have been studied in the last years in thin 
magnetic nanoelements [H [5]. The dynamical behav- 
ior of vortices has been linked to their topology [H |B]- 
[TU] . The issue of vortex dynamics has been dramati- 
cally raised anew by recent experimental observations. 
In Ref. [TT] a spin-polarized current injected into a mag- 
netic element generated a microwave signal, attributed 
to a spontaneously nucleated VA pair in rotational mo- 
tion [12] . Spin-tranfer-driven oscillations due to an ap- 
parently nonuniform magnetic state were measured in 
[131 114] . Synchronized motion of vortices could make 
the system a nanosized microwave generator [TS] . Vortex 
states have been proposed in connection to experimental 
data by micromagnetic simulations |16l 117] . In another 
set of experiments switching of the vortex polarity was 
obtained [18 21 , via dynamic generation of a vortex- 
antivortex (VA) pair [22TJ24] . 

We will show by analysis and computations that a VA 
pair where the vortex and the antivortex have opposite 
polarities, under the influence of a spin-polarized cur- 
rent, is set in a stable rotational motion due to two in- 
dependent forces: the internal interaction forces and an 
in-plane external field. The motion is stabilized by the 
spin-torque with in-plane spin-polarization. We interpret 
the obtained frequency values through a virial relation. 
The detailed theoretical study of the dynamics of nontriv- 
ial magnetization configurations is crucial for designing 
the current intensive experimental research. We further 
expect that understanding frequency generation using ro- 
bust configurations, like magnetic vortices, may lead to 



the design of devices. 

The outline of the paper is as follows. In Sec. |TT]we in- 
troduce the Landau-Lifshitz-Gilbert-Slonczewski (LLCS) 



equation. In Sec. Ill we give a rotating VA pair solution 
of the pure exchange model. In Sec. |IV| we give a virial 
relation which contains the angular frequency of rotation. 
In Sec. [V] we give numerical results for rotating VA pairs 
when spin-polarized current is flowing through a nano- 
aperture. In Sec. |VI| contains our concluding remarks. 
In the Appendix we give numerical results used in the 
figures. 



II. THE 

LANDAU-LIFSHITZ-GILBERT-SLONCZEWSKI 

EQUATION 

Our main tool is the Landau-Lifshitz-Gilbert- 
Slonczewski (LLGS) equation [35] which describes the 
dynamics of the magnetization m — (mi, 7712, ms,) in fer- 
romagnets including the effects of energy dissipation and 
spin-polarized current |26j I27| . The rationalized form of 
the equation reads 

rn = —m x {at\f — c\2Pp) — m x [m x (a 2 / + a\j3p)] , 
f = Am-m 3 e 3 + h C xt- (1) 

The dissipation constant is a and a% = 1/(1 + a 2 ), a-i = 
aa.\, the first term in f is the exchange field, in the 
second (anisotropy) term e 3 is the unit vector along 
the third magnetization direction so it describes easy- 
plane anisotropy in the (x, y) plane which is considered 
to be the plane of the film, and h cyA is an external 
magnetic field. The easy-plane anisotropy is considered 
mainly as a simplified model for the magnetostatic in- 
teraction. This is an approximation for thin films and 
it cannot be a substitute for full micromagnetic simu- 
lations. The magnetization m and the field h cxt are 
measured in units of the saturation magnetization M s , 
so we have m 2 = 1. We consider a simplified form of 
the spin-torque term where (3 is a constant, and p is 
the unit vector in the polarization direction of the cur- 
rent. We have j3 — J C /J P7 J P — noM 2 \e\df/h where 
df is the film thickness and J c is the current density. 



As an example, for a permalloy film with df — 5 nm we conjugate. The equation of motion reads 



have J p = 4.54 x e A/m . The unit of length (ex- 
change length) in Eq. (Til) is £ ox = y / 2A/(fi M^) where 
A is the exchange constant, and the unit of time is 
To = l/(^fioM s ) where 7 is the gyromagnetic ratio. For 
permalloy we have £ cx = 7.00 nm, To = 6.96 psec imply- 
ing a frequency /o = 22.9 GHz. 

Eq. (fTl) has static vortex solutions for /3 = 0, h cxt = 0. 
A vortex is a magnetization configuration characterized 
by a winding number k = 1 while an antivortex has 
K = — 1. A further topological invariant is the skyrmion 
number which is defined as the integral over the whole 
plane 



N = 



1 

4tt 



■, d x. 



-e M „(<9„m x d^m) ■ m, (2) 



where //, v — 1,2, e^ v is the totally antisymmetric ten- 
sor, and n is called the local vorticity |SJ. For a vortex 
the skyrmion number is Af — — k\/2 where A — ±1 is 
the vortex polarity. Consequently, a VA pair where the 
vortex has negative polarity and the antivortex positive 
polarity has Af = 1. Despite their nonzero Af such VA 
pairs can be generated by finite energy processes, such as 
that numerically observed in Ref. [17] in Fig. 13c. 

We will consider a thin film (or a sufficiently large mag- 
netic element) with boundary condition 



m(\r 



I -+ 00) = (1,0,0) 



(3) 



consistent with a VA pair configuration. Such a boundary 
condition may be due to the shape of the film or it may 
be imposed by the application of an external magnetic 
field of the form 



frext = (/i GX t,0, 0). 
For the spin-current term we use 

pp = (1,0,0), /3<0, 



(4) 



(5) 



which favours m = ( — 1,0,0) opposite to the magneti- 
zation at the boundary. Such polarization was used in 
Ref. [TTJ at a central region of the element and genera- 
tion of a VA pair of opposite polarities was observed. 



III. THE ISOTROPIC MODEL 

We will first investigate Eq. (fTl) for zero anisotropy, 
that is we set / = Am + h oxt . A convenient formulation 
of the problem is obtained using the complex variable 



X = 



m 2 + im 3 
1 — mi 



(6) 



where the real and imaginary parts give the stereographic 
projection of m from the point m = (1,0,0). We con- 
sider X = X{z 1 z) where z = x + iy is the complex posi- 
tion on the (x,y) plane of the film, and z is its complex 



(i-a)X = -4d z d^X+ 



-—=d z XdzX-{h^-i$)X. 

It aa 

(7) 
For /i cxt =0, j3 = any function X = X(z) is a static 
solution |28j . The simple solution X = z/ao, where ao is 
a complex constant, respects the boundary condition pi. 
It represents a vortex with negative polarity centered at 
-joo (where m 3 = —1) and an antivortex with positive 
polarity centered at iao (where TO3 = 1), and we have 
Af = 1. The constant |ao| gives vortex core size while 
the distance between the vortex and the antivortex is 

rf = 2|a | uni- 

Let us now consider nonzero and uniform ft, ext and /3. 
For X(t = 0) = z/a Eq. |7| has the solution 

*(*,«) = 7-Jr, (8) 



a(t) = a exp - oc\ [{ah cxt + /3) + i(h e 



*£)] t . 



The distance between vortex and antivortex is d(t) = 
2\a(t)\ — 2 1 ao I exp[ai(a/i ex t + P)t], while the vortices 
rotate around the origin at an angular frequency u> — 
—ct\ (ft ox t — a/3). The condition ah ext + ft — gives bal- 
ance of radial forces for, e.g., h ex t > and j3 < 0. For 
this interesting special case we have a rotating VA pair in 
steady state with an angular frequency u> — — ft ex t- For 
the case X(t = 0) = z/ao which has Af = —1 we obtain 
rotation in the opposite sense. So we conclude that the 
angular frequency is 

w = -h cxt /Af, (9) 

for Af = ±1. We have thus established that a VA pair is 
set in steady state rotational motion due to external field 
and spin-torque. Notably, the in-plane field h cxtl typi- 
cally expected to induce magnetization precession around 
its direction, is actually giving rotation of a magnetiza- 
tion configuration with Af = 1 around the axis perpen- 
dicular to the film. 



IV. A VIRIAL RELATION 

Let us now return to the full Eq. (fTl) and assume the 
existence of a rigidly rotating VA pair in steady state. 
This is expressed by the relation 



m = —u e\ v x\d v m, 



(10) 



1,2 and e\v is the totally antisymmetric 
(IT]) which gives virial 



where A, v 

tensor. Eq. (|T0| is inserted in Eq 



(integral) relations. For a uniform magnetic field Q and 
spin-torque polarization faj a, so-called, Derrick relation 
is obtained 130 



a 



^ I e + o / e ^ x \ x n d n» d x 



W a + ftext Ml + 



wiiTn CI X 



(11) 



where W a — \^{m^) 2 d 2 x is the anisotropy energy, 
dfj, v = dfj,m ■ d v m, t m = — /3(m x d^m) ■ p, and the 
integrals extend over the whole plane. The quantity 



1 = 



1 



2 j2 

3 nd x, 



P 



''2- 



(12) 



gives a measure of the size of the VA pair and is identified 
with the angular momentum 9 . We have also defined 



Mi 



x^d^m-i d x — / (1 — mi) d x, (13) 



where the last equation derives from a partial integration 
assuming vanishing boundary terms, and the last quan- 
tity gives the total magnetization (spin reversals) in the 
negative x direction. 

Let us verify Eq. (fTTl) in the case of the exchange 
model, where we should set W a — 0. For the rotating 
solution ([8]) we have J x^t^ d 2 x = and d±2 = cfei = 
0, dn =dii thus t\ v x\x^d^ v = 0. The Derrick relation 
is now greatly simplified and gives the angular frequency 
as 



u> 



We have 



4a 2 p 2 



(p 2 + a 2 ) 



2 )2 



Mi 



(2-wpdp) = pi 



(14) 



(15) 



where the first form for p\ in Eq. ( 13 ) is used. The 



integral in Eq. ( 15 ) diverges if it extends over the whole 
plane. However, the Derrick relation is valid also when 
the integrations in Eq. (11) extend over finite regions, 



except some boundary integrals are added which go to 
zero for large integration regions. We finally conclude 
that Eq. (14 1 gives u> = —h cxt for the solution pL as 
expected. 



It is useful to note that the form ( 14 ) can be obtained 



for configurations which have the following symmetries: 
mi is even in x and y, m,2 is odd in y and even in x, and 
m3 is odd in x and even in y (also, when 7712, to 3 exchange 
symmetries). We call VA pair configurations with such 
symmetries symmetric VA pairs. 



V. SPIN-CURRENT THROUGH A 
NANO- APERTURE 

We proceed to apply the ideas developed so far to an 
experimental setup where spin-polarized current is in- 
jected in a nanoelement through a nano-aperture |llj . 
Numerical simulations have shown that a VA pair of op- 
posite polarities is spontaneously created and its rotating 
motion generates microwave frequencies. The permal- 
loy nanoelements of thickness 5 nm and have an elliptic 
shape. The shape apparently favors alignment of the 
magnetization along the major axis of the ellipse effec- 
tively imposing, to a rough approximation, boundary 
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FIG. 1: VA pair in steady state rotation for spin-polarized 
current through an aperture with parameter set ( 16 1 and ap- 
plied field /i ex t = 0.05. We show snapshots at times (from 
top) t — 0, 23, 46, 69, while the period of rotation is T ~ 91.5. 
Left column shows the vector (7711,7112). Right column shows 
contour plots for 7713 = ±0.1, ±0.3, ±0.5, ±0.7, ±0.9 with pos- 
itive values represented by solid lines and negative values by 
dashed lines. We have an apparent precession of the magneti- 
zation in the vicinity of the vortex and the antivortex as they 
rotate (clockwise). 



condition similar to ([3]). The external field used is in 
the range |i? cx t| < 25 mT which, in the units of Eq. (IT]), 
reads |/i cxt | < 0.03. (since p a M s = 0.817T). A current 
/ = —4 mA through an aperture of diameter d a — 40 nm 
corresponds to a current density J G = 3.18 x 10~ 12 A/m 
and /3 = —0.70 (but partial polarization of the current 
should effectively give a smaller \j3\). 
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FIG. 2: (a) Angular frequency of rotation uj a s a function 
of the applied field h cxt for parameter set (16 1. Numerical 



results are shown by points connected by a line. The dotted 
line u) — /t cx t is plotted for comparison. Rotation is clockwise, 
(b) The vortex-antivortex separation distance d as a function 
of the applied field /i cx t for parameter set (16 1. Numerical 



results are shown by points connected by a line. 

We have simulated numerically Eq . (fl| ) using an exter- 
nal field Q and a spin-torque term (I5J which is nonzero 
only in the circular region of an aperture with diameter 
d a . We will present results for the parameter values 



13 = -0.20, 



a 



0.02, 



64 



(16) 



An initial VA pair configuration evolves under Eq. (TTj) 
and it relaxes to a steady state rotating VA pair for a 
range of values of the external field h cxt - Fig. [I] shows 
snapshots of a rotating VA pair in steady state for the 
parameter set (16) and h cxt — 0.05. We have an ap- 



parent precession of the magnetization in the vicinity of 
the vortex and the antivortex. It is important to realize 
that this process is absolutely smooth, and thus perfectly 
realizable. 

This situation should be contrasted to rotating VA 
pairs in the conservative Landau-Lifshitz (LL) equation 
(a, j3, /text = 0) which are apparently unstable states. In 
the presence of dissipation they shrink until they degener- 
ate to a point. In summary, a VA pair in the conservative 
model rotates due to the interaction between the vortex 
and the antivortex, and the rotational dynamics is linked 
to the topology of the VA pair, i.e., its nonzero skyrmion 
number |12j . 

For the full LLGS equation, solution (Isl) indicates 
that we have rotation of a VA pair due to in-plane ex- 
ternal field and the spin-torque term. This source of 
rotation is additional to the rotation due to interac- 
tions. Fig. [2k shows the angular frequency of rotation for 
< ft-cxt < 0.40. Fig. [2Jd shows the VA separation dis- 
tance d as a function of h ex t and the dependence is seen 
to be weak. As vortex positions we assume the points 
where m 3 = ±1. These are located outside the aper- 
ture for all values of h ext simulated. We do find a steady 
state also for h ext — 0, and this has a particularly large 
/ii compared to VA pairs for h cxt > 0. For h cxt > 0.40 
we do find rotating VA pairs but the radius of rotation 
fluctuates and a simple steady state is not reached [T7] , 



FIG. 3: Angular frequency of rotation u as a functio n o f 
the VA pair angular momentum i, for parameter set \1Q\ . 
Numerical results are shown by points connected by a line. 
The dotted line /xi = I is plotted for comparison. 



We actually find fluctuations also for /j cxt = 0.38, 0.40 
in some simulations (depending on the initial condition) , 
and this indicates hysteresis. 

Fig. [3k shows the angular frequency w as a function 
of the angular momentum of the VA pair I. For well 
separated vortices we expect [H] 



M^. 



(17) 



We find that formula (17) gives indeed a good approxi- 
mation to the data. 

In order to derive a guide for the understanding of the 
angular frequency u> we return to the Derrick relation 
( 11 ). The latter is satisfied with an accuracy better than 
1% in all our numerical simulations. The numerical data 
show that the term due to dissipaton (second term on 
the lhs) and the term containing f3 (third term on the 
rhs) are negligibly small in all cases. This may be due 
to the relatively small asymmetry of the simulated VA 
pairs. The Derrick relation gives for uj the approximate 
formula 



u 



i 



Mi 



(18) 



For a free VA pair (h cxt ,f3 = and a — 0) we have 
u = —W a /£ [12]. The anisotropy energy W a for a static 
isolated vortex in the LL equation (assuming no external 
fields) is equal to 7r/2 [3T], so we expect W a = 7r for a vor- 
tex pair. Using the formula (17) for the angular momen- 
tum we obtain the approximate relation u ~ —Af(2/cP) 
(for h cxt = 0). 

For understanding the angular frequency generated 
by the vortex dipole we schematically consider the two 
contributions uj — Suji + SuJ2 in Eq. (18). We find 
2.5 < W a < 3.5 for the range of /i ext simulated, so we 
argue that the contribution 5u>x = W a /£ is inversely pro- 
portional to £, as in the case of free VA pair rotation. 
The contribution 8uj2 of the second term on the rhs de- 
pends on the relation of [i\ versus I. This is shown in 
Fig. [3b, where we see that \i\ grows faster than linear 



with £. It suggests a contribution to w due to the exter- 
nal field which is 6u>2 < ft-ext for t < 145 while <5w2 > h cxt 
for £ > 145. For a different interpretation of Fig. |3ji 
consider the following. From the numerical data we see 
that 8u2 is the main contribution except for small values 
of h ex t (he., large £). One could thus take u)/h cy ± as an 
approximation for the ratio fj-i/l. 

We have performed a set of numerical simulations 
where we keep the applied field constant and vary the 
spin current strength. We briefly report here numerical 
results for 



hext = 0.05 a = 0.02, d a = 6. 



(19) 



We find steady state rotating pairs for |/3| > 0.016. For 
decreasing values of |/3| the VA distance d decreases and 
consequently the angular frequency increases. We find a 
maximum u — 0.14 and minimum VA distance d — 5 (for 
the minimum \j3\ value). Note that the d be smaller than 
d a in this case. 

Finally, we comment on the rotational dynamics in the 
presence of the magnetostatic field. We discuss briefly 
how the picture in Fig. [I] would be modified. An impor- 
tant point is that the vortices in the first and third en- 
tries (with radial orientation of the magnetization vector) 
would have higher magnetostatic energy. It is straight- 
forward to realize that the rotation of vortices would be 
accompanied by a rotation of the magnetization vector so 
that the vortex would tend to approach azimuthal mag- 
netization orientation in its vicinity. As the precession 
between radial and azimuthal vortex happens every half 
cycle, we would have the production of a frequency twice 
that of the rotating state discussed in this paper. In the 
general case, we expect that both frequencies u and 2w 
would appear in the spectrum, as measured in various 
experiments, e.g., in Refs. [LUfLl] . 

In Ref. [T7] rotating VA pairs in a full micromagnetic 
model are presented and they are accompanied by a 
satellite VA pair which is created and annihilated dur- 
ing a cycle of the main pair (the mode is termed L{). 
This process is directly compared to the comments of 
the previous paragraph. Berkov et al use numerical 
data to apply model ( 18 1 for h cxt = and find that 
it predicts a frequency 6fi — 0.44 GHz. Eq. [18] indi- 
cates that a contribution approximately 6J2 = h ex t /o 
should be added to the previous number. Using the pa- 
rameters of Ref. [T7] we have h cxt — 20mT/ (/j, M s ) = 
0.00373 and / = -yn M s /(2n) = 22.54 GHz, which 
give 5/2 = 84 MHz. The frequency is expected to be 
/ = 2(6 f 1 +6/2) = 1-05 GHz in good agreement with the 
reported value (1.2 GHz) from the simulation. 



is set in rotational motion due to two independent forces: 
the internal interaction forces between vortex and an- 
tivortex and the external field, while the motion is sta- 
bilized by the spin-torque of polarized current. Both the 
polarization of the current and the applied magnetic field 
are in-plane. We have given analytical and numerical re- 
sults for VA pairs which are in a stable steady-state rota- 
tion. We believe that our results can be used as a frame- 
work for the description of frequency generation by topo- 
logical solitons under probes, particularly, spin-polarized 
current. 

A significant step forward would be obtained if one in- 
cludes the magnetostatic field in Eq. (fTl) thus making the 
model more realistic. We expect rotational dynamics of a 
significantly more complex nature, but the gross feature 
would still be captured by the present framework. The 
magnetostatic interaction decouples rotations in physical 
space by rotations in the magnetization space and could 
thus reveal phenomena not manifested by the present 
model. 
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Appendix A: Numerical results 

We give in Tablelllpart of the numerical results used to 
plot Figs.[5][3j The angular frequency is calculated using 
the Derrick relation (TTTl) (it is in excellent agreement with 
the frequency directly measured in the simulations). 
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VI. CONCLUSIONS 

In conclusion, we have studied the dynamics of vortex- 
antivortex pairs of opposite polarities using the Landau- 
Lifshitz-Gilbcrt-Slonczewski equation [25]. The VA pair 



TABLE I: Results of numerical simulations for spin current 
through an aperture and parameter set ( 16 1 (/3 = —0.20, a = 
0.02, d a = 6). We show the distance between the vortex and 
the antivortex d, the angular momentum I and the total mag- 
netization jUl. 
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